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Abstract 

We analyze the direct CP violation in the rare decay Kl — ► 7r°e + e _ with QCD 
effects taken into account consistently in the next-to-leading order. We calculate 
the two-loop mixing between the four-quark AS" = 1 operators and the operator 
Qiv = (sd)y-A(ee)y in the NDR and HV renormalization schemes. Using the 
known two-loop anomalous dimension matrix of the four-quark operators, we 
find that the coefficient Cjy([i) depends only very weakly on [i, renormalization 
scheme and A^g. The next-to-leading QCD corrections enhance the direct CP 
violating contribution over its leading order estimate so that it remains domi- 
nant in spite of the recent decrease of |V^fe/V c b| and |U c b|. We expect typically 
BR(Kl — ► ir°e + e~)dir ~ 6 • 1CT 12 , although values as high as 1CT 11 are not yet 
excluded. 
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1 Introduction 



A clear cut observation of direct CP violation, the violation of CP symmetry in the 
decay amplitudes remains as one of the central targets of high energy physics in the 
90's and in the beginning of 21st century [pi 0, H]. The CP violation in the K — > tttt 
decays discovered almost 30 years ago can be accommodated in the standard model as 
a result of complex phases present in K°-K mixing and is usually called the indirect 
CP violation. 

The Kobayashi-Maskawa description of CP violation M predicts also the existence of 
this phenomenon in the decay amplitudes with its size strongly correlated with the weak 
couplings of the top quark and with its mass. In this respect, a very special role is played 
by the ratio e'/e. A measurement of a non-zero Re(e'/e) would automatically signal 
a direct CP violation. The experimental situation is not conclusive however. While 
the result of NA31 collaboration at CERN with Re(e'/e) = (23 ± 7) • 10" 4 clearly 
indicates direct CP violation, the value of E731 at Fermilab, Re(e' /e) = (7.4±5.9)-10~ 4 
|]] is compatible with superweak theories in which e'/e = 0. Hopefully, in about five 
years the experimental situation concerning e' / e will be clarified through the improved 
measurements by the two collaborations at the 10~ 4 level and by experiments at the 
$ factory in Frascati. 

On the theoretical side a considerable progress has been made by calculating the 
short distance Wilson coefficients of the operators contributing to e'/e beyond leading 



logarithmic approximation 0, g, || [K| |TT], [12], [13], [14]]. Unfortunately, there exist siz- 



able uncertainties in the hadronic matrix elements of these operators which hopefully 
will be reduced when new data are available. Moreover, strong cancellations between 
QCD penguin and electroweak penguin contributions for large m t make a precise theo- 
retical prediction for Re(e'/e) even harder. All efforts should be made to improve this 
situation. 

Whereas in K — > tttt decays the CP violating contribution is a tiny part of the full 
amplitude and the direct CP violation is expected to be at least by three orders of mag- 
nitude smaller than the indirect CP violation, the corresponding hierarchies are very 
different for the rare decay Kl — > TT°e + e~ . At lowest order in electroweak interactions 
(single photon, single Z-boson or double W-boson exchange), this decay takes place 
only if CP symmetry is violated J15, pi], [17], [18], [HI- Moreover, the direct CP violating 
contribution might be larger than the indirect one. The CP conserving contribution 
to the amplitude comes from a two photon exchange,[] which although higher order in 
a could be sizable. The studies of the last years |2D], [HJ indicate however that the CP 
conserving part is significantly smaller than the direct CP violating contribution. 

The size of the indirect CP violating contribution will be known once the CP con- 
serving decay Kg — > 7r°e + e~ is measured. On the other hand the direct CP violating 
contribution can be fully calculated as a function of m t , CKM parameters and the 
QCD coupling constant a s . There are practically no theoretical uncertainties related 
to hadronic matrix elements in this part, because the latter can be extracted from 
the well-measured decay K + — > tt°c + v. In what follows, we will concentrate on this 



1 Also the QED correction to a single Z-boson or a W-box exchange can give CP conserving con- 
tributions of the same order as the two photon exchange. 
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contribution relegating the discussion of the other two contributions to the end of the 
paper. 

The aim of this paper is to construct the effective Hamiltonian for Kl — > n°e + e~ , 



H. 



eff 



72 



J2 Ci{ti)Qi(ti) + C 7V (fi)Q 7V (fi) + C 7A (M W )Q 7A (M W ) 
u=i 



with the Wilson coefficients Cj(//) including leading and next-to- leading QCD correc- 
tions. Here Qi,2 denote current- current operators, Q3-6 QCD penguin operators and 

Q 7V = (sd) v _ A (ee) v , Q 7A = (sd) v _ A (ee) A (1.2) 

are the operators originating in the 7- and Z°-penguin and box diagrams of fig. |l|. For 
simplicity the CKM parameters have been suppressed in (|1.1|) . We will include them 



explicitly in section |4.1| . Our paper can be considered as a generalization of the analyses 
[ p~5| , |I6| , [T7|, [L^, [HJ to include next-to-leading logarithmic QCD effects. Whereas in 
refs. |L5|, [16], [T^, [Tj| [19| the logarithms at(a s t) n with t = ln(M^ / fi 2 ) have been summed, 
the present analysis includes also the summation of the logarithms a(a s t) n for which 
two-loop anomalous dimensions of Qi are necessary. In this way we remove the sizeable 
renormalization scheme dependence of the existing leading logarithmic results. The 
numerical significance of such a renormalization scheme dependence has been pointed 
out in ref . |2^ . 

The important observation made in refs. JT7, 0, [19| is the presence of Z°-penguin and 



box-diagram contributions which although small for mt < Mw, compete for m t > Mw 
with the photonic penguins. The contributions of these new diagrams have been already 
included in refs. [17, IS, l!|. It should be stressed however that these additional 
contributions are from the point of view of the renormalization group at the next-to- 
leading order (NLO) and their consistent inclusion requires also the summation of next- 
to- leading logarithms. In view of the NLO analyses of refs. |7j, |S], [| |It], [TTJ, |T2|, |T3], [TJJ], 
the corresponding NLO analysis for Kl — > vr°e + e _ is now substantially easier than it 

|. ^2| were performed. 



could have been at the time the calculations of refs. [17, 18 



With the complete calculation of next-to-leading contributions we will be able to 
study for the first time the dependence of the direct CP violation in Kl — > 7f°e + e~ 
on the QCD scale parameter A^jg and to remove certain ambiguities present in the 



analyses of refs. |17|, [18], |1^, ^2 



Our paper is organized as follows. In section [2], we present general formula? for 
Ci(n) beyond the leading logarithmic approximation and we discuss some aspects of 
the renormalization scheme dependence. In section ^, we calculate the two-loop mixing 
between Qi (i — 1, ... ,6) and Q 7 v, which together with the two-loop results for Qi 
of refs. [0, H, [1 HH' 0' HH HH gi ves the necessary two-loop anomalous dimension 
matrix for Kl — > 7r°e + e~ . In section |], we present the effective Hamiltonian relevant 
for Kl — > 7r°e + e _ and give numerical results for Cj(//) in the NDR and HV schemes. In 
section ||, we calculate BR(K L — > 7i°e + e~) c ii r as a function of m t and A^jg incorporating 
the updated values of the CKM parameters. We also give an analytic expression for 
BR(Kl — > Tt Q e + e~)dir using the method of Penguin-Box-Expansion |23 |. In section 
H we compare the contribution of the direct CP violation to Kl — > 7c°e + e~ with the 
other two contributions. In section a brief summary and outlook are given. A few 
technical details of the analysis have been relegated to the appendices. 
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2 Basic Formulae for Wilson Coefficients 

2.1 Operators 

Our basis of operators is given as follows 



Qi = 


\SnUa)v A[UadrJ)v A 

\" Oi. w p / V — six^fJ^Oi/ V — jtL 


Q 2 = 


(su) V - A {ud) v _ A 


Q 3 = 


(sd) v ^ A J2 q (qq)v-A 


Q 4 = 


{s a df3)v~A Sg(9/3<?a)y-A 


Q 5 = 


(sd) V - A T, q (qq)v+A 


Q 6 = 


(s a dp)v-A J2 q {qpqcx)v+A 


Q'tv = 


(a/a s )(sd) V - A (ee) v 


QlA = 


(sd) V - A (ee) A 



(2.1) 



where a and f3 denote colour indices (a, (3 = 1,...,N). We omit the colour indices 
for the colour- singlet currents. Labels (V ± A) refer to 7^(1 ± 75). The factor a/a s 
in the definition of Q' 7V allows to make all the elements of the anomalous dimension 
matrix be of the same order in a s . At the end of the renormalization group analysis, 
this factor will be put back into the Wilson coefficient C 7 y(/i) of the operator Q 7V in 
eq. (|1.2j ). There is no need to include a similar factor in Qj A as this operator does not 
mix under renormalization with the remaining operators. The sums in eq. fl2.1p run 
over the flavours active at a given scale /i. 
For /i > m c two additional current- current operators have to be taken into account 

Q\ = {s a cp) v _ A (cpd a ) v _ A , Q c 2 = {sc) V -A{cd)v-A- (2.2) 

It should be stressed that the form of Qi, Q2, Q\ and Q c 2 here differs from the one 
used by Dib, Dunietz and Gilman [0, [18[ and by Flynn and Randall where the 
corresponding Fierz conjugates have been adopted 

Qi = (sd) v ^ A (uu) v ^ Al Q 2 = {s a dp)v-A{upu a )v-A, (2.3) 

QI = (sd) v ^ A (cc) v _ A , Q c 2 = {s a dp) v „ A {c p c a ) v _ A . (2.4) 

Both bases are equally good although we somewhat prefer the basis (|2.1|), because Q2 



is taken in the colour singlet form as it appears in the tree level Hamiltonian. We will 
return later to the comparison of next-to-leading order calculations performed in these 
two bases. 

We do not include in our analysis the electroweak four-quark penguin operators 
denoted as Q7-Q10 in ref. [0 because their Wilson coefficients and matrix elements 
for Kl — > 7T°e + e~ are both of order 0(a) implying that these operators enter the 
amplitude A{K L — > n°e + e~) at 0(a 2 ). This should be distinguished from the case of 
e'/e. There, in spite of being suppressed by a/a s relative to QCD penguin operators, 
the electroweak penguin operators have to be included in the analysis because of the 
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additional enhancement factor ReA /ReA 2 — 22. Such an enhancement factor is not 
present here and the electroweak penguin operators can be safely neglected. 

Concerning the Wilson coefficients, the electroweak four-quark penguin operators 
affect through mixing under renormalization the coefficients C^-Cq at 0(a) and CVy 
at 0(a 2 ). Since the corresponding matrix elements are O(o) and 0(1) respectively, 
we again obtain a negligible 0(a 2 ) effect in A(K L — > 7r°e + e~). 

We also neglect the "magnetic moment" operators which play a crucial role in the 
b — > 57 and s — > d'y transitions. These operators, being of dimension five, do not 



influence the Wilson coefficients of the operators in eq. ( |2.1|) . Since their contributions 
to Kl — > vr°e + e~ are suppressed by an additional factor m s , they appear strictly 
speaking at higher order in chiral perturbation theory. 

2.2 Renormalization Group Equations 

The renormalization group equation for C(fi) is given by 



where (3(g) is the QCD beta function 



C(^g) = i(g)C(^g) (2.5) 



"fa) = -«-i&- A (ife — - (2 - 6) 

with 

2 38 

A) = 11 - 3/, /h = 102 - —f, (2.7) 

and / = u + d denoting the number of active flavours, u and d being the number of 
li-type and d-type flavours, respectively. 

In what follows, we will neglect the running of the electromagnetic coupling constant 
a. This is a very good approximation because only scales 1 GeV < // < Myy are 
involved in our analysis. In the numerical analysis we will take a = 1/128. For the 
effective QCD coupling constant we will use 



A)hi(Q 2 /A 



1 - 



ft In ln(Q 2 /A^) 



(2.8) 



with Af in the MS scheme. Demanding the continuity of a^'(Q) at quark thresholds 
in the form 

af\m c ) = af\m c ) } af\m h ) = af\m h ) (2.9) 

gives relations between the various Af. In what follows we will denote Aj^g = A 4 = 

^QCD- 

Although it has recently become popular to use a s (Mz) instead of Aqcd to 
parametrize the experimental knowledge about a s , we will present our results for 
the Wilson coefficients as functions of A4 in order to facilitate the comparison with 
refs. |T7], [TH], [TjJ. For each particular value of A 4 used, we give the corresponding value 
of a s (Mz) with the next-to-leading order accuracy in app. |D|. 



4 



Since the operator Qja does not mix with the remaining operators in eq. ( |2.1| ), we will 
not include its coefficient in C(\i). Consequently, ^(g) in eq. ( |2.5|) is a 7 x 7 anomalous 
dimension matrix which we expand as follows 



nyj 7 16vr 2 y (16tt 2 ) 2 
The calculation of this matrix will be presented in section 



2.3 Solution for C(fi) 

Following we write the solution of eq. ( |2.5| ) as 

C(fi) = U(p, M W )C(M W ) 

where the renormalization group evolution matrix is given generally by 

rg(mi) A T (g') 
U (mi, ma) = T g exp / dg 

J g ( m2 ) (3(g') 



(2.10) 



(2.11) 



(2.12) 



with 777,1 < m 2 . Here T g denotes such ordering in the coupling constants that they 
increase from right to left. Expanding j(g) and (3(g) as in eqs. ( |2.1(J| ) and (|2.6|) respec- 
tively one finds 



U(m^m 2 ) = (i + 2Mj U^( mi ,m 2 ) (l - 



(2.13) 



where tA°) (777,1, 777,2) denotes the evolution matrix in the leading logarithmic approxi- 
mation and J summarizes the next-to-leading corrections to this evolution. If 



7§> ee V^ 7 (0)T t> 



G = V-^ T V 



(2.14) 



where 7^ denotes a diagonal matrix whose diagonal elements are the components of 
the vector 7^, then 

U (0) (m u m 2 ) = V 
For the matrix J one gets 



' a s (m 2 ) \' 
a s (mi) ) 



V 1 with a 



7 



(o) 



D 



2[3 



(2.15) 



J = vsv- 1 

where the elements of S are given by 

(0) A 



Sij — <%7j 



2/5o 2 2A, + 7f } - 7 7 (0) 



(2.16) 



(2.17) 



with 7 i - ' 1 denoting the elements of 7^ and Gij the elements of G. 
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2.4 The Initial Conditions C(M W ) 

The calculation of the initial conditions C(Mw) has been discussed at length in refs. |L7], 
18, 19, 12, 13 . In fact, the results for C(Mw) can be extracted from these papers. 
Therefore we recall here only the basic ingredients of this calculation stressing one 
point which has not been discussed in previous analyses of K L — > ir°e + e~ . 

In order to find C(M\v), the one-loop current-current and penguin diagrams of fig. [3] 
with the full W and Z propagators and internal top quark exchanges have to be calcu- 
lated first. Also the relevant counterterms have to be added to make the result finite. 
Subsequently, the result of this calculation has to be expressed in terms of matrix el- 
ements (Q(M W )). The latter are found by inserting the operators Q as depicted in 
fig. [| into the one-loop current- current and penguin diagrams of fig. || and calculating 
the finite contributions in some renormalization scheme. 

The resulting coefficients C(Mw) do not depend on the assumptions made about 
the properties of the external lines in figs. [I] and [3], i.e. on the infrared structure of 
the theory. They depend however on the renormalization scheme through the matrix 
elements (Q(Mw))- In our calculation we will use as in ||, |9], [U], [IT], [L2|, |IJ| [IJJ two 
renormalization schemes usually called NDR and HV. 

For the NDR scheme and the basis ( |2.1| ) one finds 



Ci(M] 



w) 



H« S (M, 



w i 



C 3 (M, 



w / 



47T 



-C 4 {M W ) = C 5 {M W ) 



C 2 {M W ) = 1 - 



11 a s (M w ) 



6 4vr 



-MM, 



w ) 



a s (M, 



H ' J 



24vr 



E(x t ), 



C' 7V (M W ) 



a s {M } 



w , 



2tt 



C(x t ) - B(x t ) 



sin 2 6 



D{x t )-AC{x t ) 



vi- 



and 



Here 



C 7A {M W ) 



E(x t ) = E{x t ) - -, 



a B(x t ) - C(x t ) 
2n sin 2 9w 



D{x t ) = D{x t ) - -, 



and 



B(x) 
C{x) 
D(x) 
E(x) 



1 

4 

x 



x 



1 — x 
x — 6 



+ 



+ 



xhix 
x-1) 2 
3x + 2 



m t 



x — 1 [x — 
-19x 3 + 25x 2 

36(x-l) 3 + 
x(18 — llx — x 2 ] 
12(1 -x) 3 



I) 2 

,.2 



\nx 



x (5x — 2x — 6) 



+ 



lnx - 

18(x - l) 4 

x 2 (15- 16x + 4x 2 ) 



6(1 



x 



Unx, 



lux In x. 

3 



2.18) 

2.19) 
2.20) 

2.21) 

2.22) 
2.23) 

2.24) 
2.25) 
2.26) 
2.27) 
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Comparing these results with the initial conditions used in refs. []I7|, [TH], [T!| we notice 



some differences. In refs. [17, 18, 19] the leading order evolution matrices have been 



used but in the initial conditions the mj-dependent next-to-leading terms have been 
already taken into account. These mj-dependent terms represented by the functions 
B(xt), C(xt), D(xt) and E(xt) are renormalization scheme independent. Yet, the 
initial conditions may also contain scheme dependent next-to-leading terms which of 
course were not present in the analyses of refs. ^L7[ |18| , |l9|j . In the NDR scheme these 
scheme dependent terms enter the initial conditions for Qi and Q2 and also modify the 
functions E(xt) and D(xt) by constant terms as given in eq. (|2.22|) . 
The additional constant terms in eq. ( |2.22| ) are absent in the case of the basis ( |2.3| )- 
used in refs. |L7|, [18], |19| . However, as pointed out in ref. || |I2|1 , in this case 



also the corresponding two-loop anomalous dimensions 7^ change so that the physical 



amplitudes remain independent of whether the basis ( |2.1|) or (2.3)- (|2.4|) is used. This 
clarifies the observation made by Flynn and Randall fl9j that the initial conditions for 
penguin operators do depend on the form of the operators. This dependence signals 
the necessity of going beyond the leading logarithmic approximation for [/(mi,m 2 ) in 
order to obtain physical results. 

In the HV scheme the constant terms in eq. (|2.22|) are also absent. In addition 
the initial conditions C\^{My/) change. They can be found in app. |A]. Again, in 
physical quantities this change is cancelled by the corresponding change in the two- 
loop anomalous dimension matrix 7W calculated in the HV scheme. 



2.5 Renormalization Scheme Dependence of 

The renormalization scheme dependence of 7W and the cancellation of scheme depen- 
dences in the final physical amplitudes have been discussed at length in refs. || [| [H], 
11, 12, 13, 14|. Here we would like only to recall one important relation between (j^) a 
and ( , y^)b calculated in two different renormalization schemes a and b or two different 
bases like Q and (p|)-(^l. We have 



( 7 (1) ) 6 = (7 w )a+ Af,7 W +2p Af 



with Af 



(r) a , where (f)j are defined by 

a s (M w ] 



(Q(M w ))i 



1 + 



47T 



(0.- 



(Q (0) ). 



(2.28) 



(2.29) 



Here, is a tree- level matrix element and {Q)i denotes renormalized one- loop ma- 
trix elements calculated in the scheme i. The matrices (f)j are obtained by calculating 
the finite terms in the one- loop diagrams of fig. [3|. 

Relation (|2.28| ) is very useful as it allows to test compatibility of the two-loop anoma- 
lous dimension matrices calculated in two different renormalization or regularization 
schemes. It also plays a central role in demonstrating the scheme independence of phys- 
ical quantities || . Moreover, as demonstrated in refs. [|H], [TT| and in the app. |C] here, 
the relation ( |2.28| ) allows to avoid explicit calculation of two-loop diagrams involving 
traces Tr{pf^y v ^p^xl5) which are dangerous in the NDR scheme. 
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3 Anomalous Dimension Matrices 



3.1 The Matrix <y(°) 

The 6x6 submatrix of ■y^ involving the operators Q\-Qe has been calculated in 
ref. |L5, |24], [26|, [27] . In our normalizations it is given in app. A of ref. [10|. Here we 
only give the remaining non- vanishing entries of j(°>: 



iff = 


9 




(0) 
727 = 


16 

~Y 


iff = 


16 - 


' d 1\ 
C-2-n) 


(0) 
747 = 


16 
~~9~ 


(0) 
757 = 


9 ' 




(0) 
767 = 


16 


(0) 
777 = 


-2/3„ = 


22 4 

iV+ -/ 

3 3 17 







-N 




(3.1) 
(3.2) 

(3.3) 
(3.4) 



where N is the number of colours. These elements have been first calculated in [plj 
except that 737'* and 747 have been corrected in refs. [19, |28| . We agree with these final 
results. 



3.2 The Matrix 7 



(i) 



The 6x6 submatrix of 7W involving the operators Q\-Qe has been calculated in 
ref. R HM in the NDR and HV renormalization schemes. The results presented there 



include also the ones obtained in refs. [[?], § where the two-loop mixing of Qi and Q 2 
operators has been considered. In our normalizations, the submatrix in question is 
given in app. B of ref. [10 . 

In order to calculate the remaining entries of 7' 1 ', one has to evaluate the two- loop 
penguin diagrams of fig. |J Since these diagrams constitute a subset of the diagrams 
considered in [Kj], all the singularities necessary for the evaluation of the entries in 



question can be extracted for instance from tab. 2 of ref. [I0]|. Consequently, it was 
enough to change the relevant colour and electric charge factors in the diagrams of 
fig-i 

Our final results for the seventh column of 7^ are given in the NDR scheme as 
follows: 2 



41 



iff 



7S 1 



iff 




7§» 



81 V N 



280 ,64\ / 1 

u - + d Ti) far 




N 



N' 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



2 The first six entries in the seventh row of 7^ 1 - ) vanish. 
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iff 



N 



(3.10) 



where Cp — N 2 m X ■ The corresponding results for the HV scheme can either be obtained 
by direct calculation or by using the relation ( |2.28 ). They can be found in app. |A[ 
As shown in app. [FJof the present paper, the results (ft.5| ) -( jO0|) are in agreement with 



those calculated independently in ref . |29f , as long as the differences in the conventions 
used are properly taken into account. However, the present paper is the first one where 
the the two-loop mixing among Q\-Q% from ref. |§ is consistently put together with 
the mixing of Q\-Qe into Qiy The author of ref. [EUJ did not take into account that 
certain conventions used in ref. || were different from his own, and did not transform 
the results of ref. into his own conventions (see app. |B| for details). 

The calculation of two-loop mixings made in the NDR scheme requires some care 
when Tr(7 M 7„7 p 7 A 7 5 ) appears in the calculation. One of the ways dealing with this 
problem is just to avoid appearance of such traces by making use of the fact that a Fierz 
rearrangement of any of the considered operators should not affect physical results. 
This method has been discussed at length in refs. 



|T_TJ . It is straightforward to 
extend it to the case of the operator basis containing also Qjv- For the more complex 



cases of 737 and 747 we describe this extension in app. 0. This is the method we have 
used to obtain the results ( |3.6|) and (|3 . 7| ) . 

Another method described in app. A of ref. [21] ] amounts to treating fermion lines with 
single {7^,75} as evanescent structures, i.e. structures that vanish in 4 dimensions but 
have to be treated as independent basic structures in D dimensions. In consequence, 
no algebraic inconsistencies with anticommuting 75 arise.j] Getting rid of the trace 
7V(7/i7t/7p7A75) is achieved by introducing other evanescent structures. In the particu- 
lar calculation of the two-loop mixing among the operators of the basis ( |2.1| ), the new 
evanescent structures can be shown to give no contribution. So the calculation in our 
particular case effectively amounts to anticommuting 75 whenever necessary (except 
for traces with single 75) and using eq. (19) of ref. |29) (or, equivalently, eq. (5.20) of 
ref. pi| ) to get rid of the traces with 75. 

Actually, in refs. [10], [Ll| and in ref. the two approaches were used and checked 



to give the same results. Also the equivalence of the physical results to the ones 
obtained within the HV scheme has been checked in refs. 
of mixing among the four-quark operators. 



TO, HI H3L HI for the case 



4 The Effective Hamiltonian for Kl —* tt e + e 
Beyond Leading Logarithms 

4.1 Master Formulae for Wilson Coefficient Functions 

Integrating out simultaneously W, Z and t we construct the effective Hamiltonian 
for AS* = 1 transitions relevant for — > 7r°e + e~ with the operators normalized at 

3 The triangle anomaly is also correctly reproduced. 
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fi = M\y. Following the procedure of section [2.4| we find 



where 



H eff (AS = 1) = ^V ud V: s {(1 - r)H u \M w ) + rH ut {M w )} (4.1) 



#« C (M W ) = C 1 (M w ,)(Qi-g i ) + C , 2(M w )(g 2 -Q^), (4.2) 

# ui (M w ) = ^C J (M w )Q, + CV(M w )Q / 7y + C 7A (M w ,)g M (4.3) 
i=i 



and 



(4.4) 



Here Qi and 2 are the operators defined in eqs. (|2.1|)-( |2~2"D with q = u,d, s, c, b. Due 
to the GIM mechanism, there are no contributions of Qi{i ^ 1,2) to H uc (My/). The 
coefficients Ci(M w ), C' 7V (M W ) and C 7A (M W ) are given in eqs. (gfg)-(ggT]). 



We next integrate out successively the 6- and c-quarks and transform H e f f of eq. (|4.1|) 
from \i = My/ down to /J < m c relevant for Kl — > 7i°e + e~ . This transformation can 
be compactly described as follows. 

H UC (M W ) — E^W + ^Wtv, ( 4 -5) 

i=l 

# ut (M w ) — E^(/ u )^ + ^(^)Q7V + C , 7a(M w )Q 7A (4.6) 



with the operators and Q 7 y renormalized at the scale /x. There is no renormaliza- 
tion of Qta- Since the charm quark has been integrated out, the penguin operators 
contribute also in eq. ( |4.5|) . 

In writing eqs. ( |4.5|) and (|4.6| ) we have neglected the internal u- and c-quark con- 
tributions to the coefficient of the operator Q 7 a because due to the vanishing of C(x) 
and B(x) as i lni for x — > they are very small compared to the top contribution. 

The coefficients (zi,z' 7V ) and (vi,v' 7V ) are components of the seven-dimensional col- 
umn vectors 

z(n) = U 3 (n,m c )z(m c ) (4.7) 

and 

v{fi) = U 3 {fi, m c )M(m c )U 4 (m c , m b )M(m b )U 5 (m b , M W )C{M W ) (4.8) 

with the components of C(Mw) given in eqs. ( |2.18| )-( |2.20| ). 

The evolution matrices Uf are given by eq. ([2.13D with / denoting the number of 
active flavours. M(m») is the matching matrix at quark threshold m^. It will be 



discussed separately in section |4.2 



Inserting eqs. ( |4.5| ) and ( |4.6| ) into eq. ( |4.1| ) and introducing the coefficients 



a a 
a s {fx) a s {ii) 
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we find the effective Hamiltonian relevant for the decay Ki — > it e + e 



G F 

H eff = -^V ud V* 



TV 

Y^{zi{n) + ryi{ij))Qi{fj) + ry 7A (M w )Q 7A 



(4.10) 



where 



y;0) = Vii/i) - z^), V7a{M w ) = C 7A (M W ). (4.11) 

In order to find z(fi) we still need z(m c ). Due to the GIM mechanism, the coefficients 
Zi(ii) of the penguin operators Qi, i ^ 1,2 are zero in 5- and 4-flavour theories. The 
evolution in this range of \i {p, > m c ) involves then only current-current operators 
Q12 = Qi 2 and Q\ 2 with the initial conditions 

z x {M w ) = C X {M W ), z 2 (M w ) = C 2 (M W ) (4.12) 



where C\^{Mw) are given in eq. ( |2.18|) . Qi 2 an d Q\ 2 do not mix under renormalization, 



and the coefficients of these two pairs of operators are equal for \i > m c . We then find 

( Sk) ) = °' {m - A « ( "^tl ) (4 - 13) 

where this time the evolution matrices £4.5 contain only the 2x2 anomalous dimen- 
sion matrices describing the mixing between current-current operators. The matching 
matrix M(mb) is in this case a 2 x 2 unit matrix. When the charm quark is integrated 
out, the operators QI 2 disappear from the effective Hamiltonian and the coefficients 
Zi(n), i 7^ 1, 2 for penguin operators become non-zero. In order to calculate Zi(m c ) for 
penguin operators, a proper matching between effective 4- and 3-quark theories has to 
be made. This matching has been already discussed in detail in section 4 of ref. [|T2] 



where the first six components of z(m c ) are given. In order to find z' 7V (m c ) which 
results from the diagrams of fig. |||, we simply have to rescale the result for the z 7 (m c ) 
in eq. (4.29) of ref. |12[ by a factor of — 3a s /a. 



In summary, the components of z(m c ) are given by eq. (|4.13| ) 



and 



where 



/ \ Zi{m c ) z e (m c ) a s {m c ) 
Z3{m c ) = — = z 5 (m c ) = — = ^^ F *( mc ) ( 4 - 14 ) 

= -^4r^F e (m c ) (4.15) 



F,(m c ) = { -^ {m °>' ™ fl (4.16) 
F,W = J -|(3 2 ,K)+^K)), NDR (41?) 



In the NDR scheme these results are given for the basis (|2Tl| 
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4.2 The Matching Matrix M 

The additional matching matrix M between the operators in an /-flavour effective 
theory and an (/ + l)-flavour effective theory is necessary because the QCD penguin 
operators have different structure in different effective theories. In our case 

M(m) = 1 + -^^5fJ (4.18) 

where 8rT is a 7 x 7 matrix. 



The procedure for finding M(m) has been outlined in ref. [12] where the 6x6 
submatrix involving Q\—Qq has been calculated. The remaining entries of Sf s can be 
found from the matrix 8r e , also calculated in ref. | l2f , by making a simple rescaling by 
— 3a s /a as in the case of z(m c ). 

In summary, 5r s is given as follows. The non-vanishing rows in the 6x6 submatrix 
are given for /i = m& and ji = m c by 

(Sr s h = (Sr s h = —Pi (4.19) 
y 

where 

P= (0,0,-1,-1). (4.20) 

The seventh row of 8r s vanishes. The seventh column has the following non-vanishing 
elements for fi = nib- 

20 

(5r s ) 37 = 3(5f s ) 47 = (5f s ) 57 = 3(5f s ) 67 = --. (4.21) 
For /i = m c , — y has to be replaced by y. 

4.3 Numerical Results for Wilson Coefficients 

In tab. ^ we give numerical results for the coefficients z\ and with i — 1, ... ,6, 7V, 7 A 
in the NDR scheme, for A^g = 200, 300, 400 MeV, fj, = 1 GeV and m t = 150 GeV. 
The results for the coefficients of Q±-Qq are given for completeness. They have been 



analyzed in ||12|| . The numerical results for them presented here differ slightly from the 
ones in |12| because in the present paper the effect of electroweak penguins, being of 
higher order in a, has not been included. In tab. we show the /i-dependence of z-jy ja 
for different values of A^g and for m t = 150 GeV. The coefficient z 7V does not depend 
on m t . In tab. [| we show the //-dependence of y 7 y /a for different values of A^ and for 
m t = 150 GeV. In tab. f| we show the m t -dependence of y 7 v/ a an d ViaI^ for different 
values of A^g and for /x = 1.0 GeV. Concentrating on the coefficients of Q 7 y and Q 7 a 
we observe: 

• The coefficient y 7 y is for m t = 150 GeV enhanced by 30-40% over the leading 
order result so that the suppression of this coefficient through QCD, found in 
T7\ |H| |T9j] , is considerably weakened. Without QCD corrections one would find 
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y-jv/ a — 0.795. For higher values of m t the effect of next-to- leading order cor- 
rections, albeit sizeable, is weaker: 25-35% enhancement for m t = 200 GeV and 
smaller for higher mt- This feature will be clearly understood in the presentation 
in section 5. 

• The scheme- and /i-dependences of y 7 y are very weak as shown in particular in 
tab. [| Also the dependence of yrv on is rather weak, at most of 3% in 
the range of parameters considered. The weak /^-dependence is related to the 
fact that Qjv does not carry any anomalous dimension and that for fi < m c the 
coefficient y 7 y is only influenced by QCD penguin operators because yi = y 2 = 0. 

• The scheme- and /^-dependences of z 7 y shown in tab. are very strong. This 
implies that the calculation of the real amplitudes such as Kg — > 7r°e + e~ requires 
the inclusion of the matrix elements of Q\-Q% in order to cancel the scheme- and 
/t-dependences present in z 7 y. Since the matrix elements (7i e + e~\Qi\K s ) can 
only be calculated using non-perturbative methods, it is clear that Ks — > n°e + e~ 
and consequently the indirect CP violation in Kl — > 7r°e + e~ are much harder to 
calculate than BR{K L — > TT e + e~) c iir ■ 

• Inspecting tab. ^ and fig. ||we observe in accordance with the analyses of Jl7|, |TJ, 
pTOfl that at a certain value of m t the coefficient \yrA\ becomes larger than \yrv\- 
For LO this transition takes place roughly for m t = 140 ± 10 GeV, where the 
precise value depends on Aj^g. The next-to-leading order corrections shift this 
transition point to m t = 180 ± 5 GeV. 



5 BR(Kl —> n e + e Beyond Leading Logarithms 



5.1 Analytic Formula 

Let us introduce the numerical constant 



« = T/j^y BR{K+ ^ e+v) = 416 • (5 - 1} 



We find then as in [18 



BR(K L -> 7T°e + e-) dir = K(lm\ t ) 2 (y 2 7A + y 2 7V ) (5.2) 

where 

ImAt = Im(V td V*) w \V ub \ \V cb \ sin 5 (5.3) 

in the standard parametrization of the CKM matrix. 

In deriving eq. (|5.2|) , the matrix elements of the operators Q^-Q% have been neglected, 
i. e. we have assumed that 

^(/^(vrVe-IQ^) « y 7V ^)(7c°e + e-\Q 7V \K L ). (5.4) 

i=3 
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This assumption is supported by the corresponding relation for the quark-level matrix 
elements 

6 

^2vi(ji)(de + e~\Qi\s) < y 7V (/j,)(de + e'\Q 7V \s) (5.5) 

i=3 

that can be easily verified perturbatively. The neglect of the matrix elements of the 
QCD penguin oerators is compatible with the very weak /^-dependence of the retained 
contribution. 



Using next the method of the penguin-box expansion |23] we can write 



Vrv = n a 2 n ( P o + Y ( x t) ~ 4 sin 2 9 W Z(x t ) + P E E(x t )) , (5.6) 
2tt sin W v ' 



Via = - Y(x t ) (5.7) 

2tt sin Uy/ 



where the gauge invariant combinations of B(x), C(x) and D(x) are given as follows: 
Y{x)=C(x)-B(x), Z(x) = C(x) + ~D(x). (5.8) 



To a very good approximation one has for 120 GeV < m t < 300 GeV 

Y(x t ) = 0.315x°- 78 , Z{x t ) = 0.175x°- 93 . (5.9) 

The constants Pq and Pe summarize the contributions from scales below fi = M^y. 
Pe is 0(1O~ 2 ) and consequently the last term in eq. ( |5.6|) can be neglected. Pq is given 
for different values of fi and Ay^ in tab. |5]. We show there also the leading order results 
and the case without QCD corrections. The analytic expressions in eqs. ( |5.6|) and (|5.7p 
are useful as they show not only the explicit m r dependence, but also isolate the impact 
of leading and next-to-leading QCD effects. These effects modify only the constants P 
and Pe- As anticipated from the results of section |], Pq is strongly enhanced relatively 
to the LO result. This enhancement amounts roughly to a factor of 1.6 ± 0.1. Part of 
this enhancement is however due to the fact that for Alo — A-^ the QCD coupling 
constant in the leading order is 20-30% larger than its next-to-leading order value. 
Calculating Pq in LO but with the full a s of (|2.8|) we have found that the enhancement 
amounts then to a factor of 1.33 ± 0.06. In any case the inclusion of NLO QCD effects 
and a meaningful use of Aj^g show that the next-to-leading order effects weaken the 
QCD suppression of y 7 y- As seen in tab. [|, the suppression of P by QCD corrections 
amounts to about 15% in the complete next-to-leading calculation. 

We next find that 

f 0.332 for m t = 150 GeV 
Y(x t ) -4sin 2 ^ Z(x t ) = I 0.431 for m t = 200 GeV (5.10) 

[ 0.523 for m t = 250 GeV 

whereas in the same range of m t the function Y(x t ) changes from 0.840 to 1.863. 
Because of the additional constant term P in eq. (|5.6|) , the m t -dependence of y 7 y is 
rather weak and consequently the dominant m r dependence of BR(Kl 7r e + e~)dir 
for fixed lm\ t originates from the coefficient of the operator Q 7 a- In fact, for m t > 



14 



lSOGeV^ this operator gives larger contribution to the branching ratio than Qjv This 
is shown in fig. || where \y 7 y/a\ 2 and \y 7 A/ct\ 2 are plotted as functions of m t together 
with the leading order result and the case without QCD corrections. Finally, we note 
that for m t < 200 GeV the constant term proportional to Pq constitutes at least 60% 
of y 7V . 



5.2 Uncertainties Due to the Definition of m t 

At the level of accuracy at which we work we cannot address the question of the 
definition of m t used here. In order to be able to analyze this question, one would have 
to calculate perturbative QCD corrections to the functions Y(xt) and Z(xt) and include 
also an additional order in the renormalization group improved perturbative calculation 
of Pq. The latter would require evaluation of three- loop anomalous dimension matrices, 
which in the near future nobody will attempt. In any case, we expect only a small 
correction to Pq. The uncertainty due to the choice of /i in m t (fi) can be substantial, as 
stressed in refs. 0, [H] ^2|, and may result in 20-30% uncertainties in the branching 
ratios. It can only be reduced if 0(a s ) corrections to Y(xt) and Z(xt) are included. 
For K + — > 7T + b>v, Kl — ► 7T°vv, B — ► and B — > X s vv this has been done 

in refs. pO] [31] , |32|. The inclusion of these corrections reduces the uncertainty in 
the corresponding branching ratios to a few percent. Fortunately, the result for the 



corrected function Y\xA given in refs. [30, pi], 32 can be directly used here. The 



message of refs. |30], [H], |32[ is the following: For m t = m t (m t ) : the QCD corrections to 
Y(x t ) are below 2%. Corresponding corrections to Z(x) are not known. Fortunately, 
the m t -dependence of y 7 y is much weaker and the uncertainty due to the choice of fi in 
m t (n) is small. On the basis of these arguments and the result of refs. |3(| [31], |32| we 
believe that if m t = m t (m t ) is chosen, the additional QCD corrections to BR(K^ — > 
Tr e + e~)dir should be small. 



5.3 Numerical Analysis 



In the numerical calculation we need the value of lm\ t which can be extracted from 
the usual analysis of the parameter e that measures indirect CP violation. The formula 



for e can be found for instance in section 7 of ref. |12| as well as in many other papers 
and will not be repeated here. We only give the values of the parameters which we use 
here. We take \V US \ = 0.22, 



T] 2 = 0.570, r] 3 = 0.36|§, H, |37|, [3§ (leading order). (5.11) 



and 

Range a: |V c6 | = 0.040 ± 0.005, \V ub /V cb \ = 0.08 ± 0.02, B K = 0.70 ± 0.20, (5.12) 

The analysis of e gives two solutions for the CKM phase 5. The values of lm\ t 
decrease with increasing m t and this decrease is stronger for 5 in the second quadrant. 
The dependence of lm\ t on m t compensates the corresponding m r dependence of y 7 y 
and y-jA- As a result, BR(Kl — > Tc°e + e~)dir decreases with m t for the solution II and 
is increasing with mt for the solution I substantially slower than one would expect 
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from eq. ( [5~2"| ) would Im\ t be independent of m t . In tab. |G] we show the values of 
BR{K L -> vr°e + e-) dir in units of 1(T 12 for B K = 0.7, \V cb \ = 0.041, \V ub /V cb \ = 0.10 as 
a function of m*. For m t < 140 GeV, e cannot be fitted for these choice of parameters. 
In fig. |] we present the results when the parameters are varied in the full range given 
in (|5.12 ). In order to illustrate the impact of possible future improvements in the 
determination of B^ and the CKM parameters, in figs. [7] and |] we show the results 
when smaller ranges of parameters given by 



Range b: \V cb \ 
Range c: \V cb \ 



0.040 ±0.002, 
0.043 ± 0.002, 



IWV^I = 0.08 ± 0.01, B K 
\VjV cb \ = 0.09 ± 0.01, B K 



0.75 ±0.05 (5.13) 
0.55 ±0.05 (5.14) 



are used. B^ in range b is in the ball park of most recent lattice (Bk = 0.825 ± 
0.027 ±0.023) |3 and 1/N (B K = 0.7 ± 0.1) g(J results. In range c it is closer to the 
values obtained in the hadron duality approach (Bk = 0.4 ± 0.1) j|T|. In the latter 
case we had to take larger values for \V cb \ in order to obtain solutions in the e-analysis. 



Finally in fig. |9| we show the ratio BRiK^ 



7r°e + e 



)/(ImA 4 ) 2 as a function of m t . 



All these results are self explanatory and we only make a few comments: 

• The NLO corrections enhance the direct CP violating contributions with respect 
to its LO estimate. 

• Due to large uncertainties in ImAt which enters the branching ratio quadratically, 
the enhancement found here cannot be appreciated yet. However as seen in 
figs. [F] and |8] the more accurate data for lm\ t in the future will allow to feel this 
enhancement clearly. 



For central values of the parameters in range a and mt = 170 ± 30 GeV we find 

uu(zr o + -n / (4.9 ±0.5) ■ 10" 12 , Solution I 

BR(K L - 7r°e + e ) dir = n _ 12 ' Q , + . TT (5.15) 



but as seen in fig. |^ values as high as 10 11 are not excluded when the solution I 
is chosen. 

The analysis of B°-B° mixing can constrain further the predicted values of BR(Kl — > 
7r°e + e~ )dir- The usual box diagrams with top quark exchanges give for the mixing 
parameter 



where 



We use 



x d = 3.89 ■ 10 



S{x) 




V td \ 2 S{x t ) 




(1-x) 2 ; 2(x-l) 



'B B F B = 200 ± 30 MeV 



\nx. 



(5.16) 

(5.17) 
(5.18) 



which is in the ball park of various lattice and QCD sum rule estimates P2[ . Setting 
tb = 1.5 ps we require 

x d = 0.72 ± 0.08 [43] (5.19) 
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in accordance with the most recent average of CLEO and ARGUS data. The impact of 
B°-B mixing with \[B~bFb as given above is shown in figs. |6|-[| for the NLO results as 
a solid line within the dark shaded area. With this additional requirement the region 
to the right of the solid line is excluded. The constraint is invisible for solution I. For 
solution II we observe that for higher values of m t the constraint coming from B°-Tf 
mixing is rather effective. These features are in accordance with general expectations 
that for high values of m t and Fb the solution I is favoured and that with more precise 
high values of Fb the solution II could even be excluded. This would also constrain 
the values for BR(Kl — > 7r°e + e~ for solution I. 



6 Comparison with the Other Two Contributions 

Now we want to compare the results obtained in the previous section with the estimates 
made for the indirect CP-violating contribution and the CP-conserving one. The most 



recent summary of the present status has been given by Pich [^T| . We have nothing to 
add to this issue here and only state his final conclusions. 

The indirect CP violating contribution is given by the K$ — > 7r°e + e~ amplitude times 
the CP parameter e. Once BR(Ks — > ir°e + e~) has been accurately measured, it will 
be possible to calculate this contribution precisely. Using chiral perturbation theory it 
is however possible to get an estimate by relating K$ —>■ TT°e + e~ to the K + — ► 7i + e + e~ 
transition. Present data implies 

BR(K L -> n°e + e~) mdir < 1.6 • 1(T 12 , (6.1) 

i. e. a branching ratio more than a factor of 2 below the direct CP violating contribution 
for the CKM phase S in the second quadrant. More importantly, if the solution I is 
chosen, the direct CP violating contribution is by a factor of 5 or more larger than the 
indirect one. This should be contrasted with K — > tttt decays. It should be stressed, 
however, that in reality the CP indirect amplitude may interfere with the vector part 
of the CP direct amplitude []. This interference could have both signs. For a very heavy 
top quark for which the first quadrant is favoured and the axial-vector contribution 
to the direct CP violation dominates, the presence of the indirect CP violation will 
not have a large impact on our analysis given above. Yet for m t = 0(150 GeV) and 
5 in the second quadrant, the indirect CP violation cannot be neglected and a more 
careful estimate of the full CP violating branching ratio requires a better analysis which 
includes the interference in question. 

The estimate of the CP conserving contribution is more difficult. We refer the reader 
to ref. [EI] where the literature on this subject can be found. The most recent analysis 



in the framework of chiral perturbation theory gives [ZO] 

BR{K L -> n°e + e-) cons « (0.3 - 1.8) • 10~ 12 , (6.2) 

i. e. well below the direct CP violating contribution. An improved estimate of this 
component of Kl — > 7t°e + e~ is certainly desirable. 



4 The axial part of the CP indirect amplitude can be safely neglected 
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7 Summary and Outlook 



We have calculated the direct CP violating contribution to Kl — > n°e + e~ beyond the 
leading logarithmic approximation. This required the evaluation of two-loop mixing 
between the four-quark operators Q±-Qe and the semileptonic operator Qyv- Our main 
findings can be summarized as follows: 

• The next-to-leading QCD corrections enhance the branching ratio BR(Kl — > 
TT°e + e~)dir so that its suppression found in the leading order analysis is consid- 
erably weakened. 

• The net QCD suppression of this contribution for m t > 150 GeV is at most 20%. 

• The very weak dependences of the resulting BR(K L — > TT°e + e~) dir on and 
fj, indicate that this branching ratio can be reliably calculated in the Standard 
Model, provided m t and the CKM parameters have been better determined. 



With the help of the two-loop calculation of refs. fl30| , |3T| , |32| we have also some 
control over the renormalization scale of m t . Consequently when m t = m t (m t ) 
is chosen we expect on the basis of refs. PU|, [5T], |3"2"|] that still higher order QCD 



corrections are very small. 

Provided the phase 5 is in the first quadrant and m t = 170 ± 30 GeV we expect 

(6.0 ± 4.0) • 10" 12 , Range a 
BR{K L -> n e + e-) dir = { (4.5 ± 1.5) • 10~ 12 , Range b (7.1) 

(7.5 ± 2.5) • 10- 12 , Range c 



for the three ranges of parameters considered. As seen from (7.1) the present 
prediction, range a, will be improved when the CKM parameters and Bk have 
been better determined and the top quark has been discovered. For 5 in the 
second quadrant the branching ratio is typically by a factor of two smaller. 

These results confirm the expectations made previously by other authors that the 
decay Kl — > 7r°e + e~ is dominated by direct CP violation (see ^TJ and references 
therein). The present experimental bounds 



Bfl(AWVe-) < " S (7-2) 




are still by three orders of magnitude away from the theoretical expectations in 
the Standard Model. Yet the prospects of getting the required sensitivity of order 
10~ n -10~ 12 in five years are encouraging £| [|. 
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Appendices 



A Calculation in the HV Scheme 



The Wilson coefficients in the HV scheme can be found using 



47T 



(a.i; 



where Af has been defined in ( |2.29| ). The 6x6 submatrix of Af relating HV and NDR 
schemes is given in eqs. (3.9)-(3.13) of ref. f9j. Performing an analogous calculation in 
the presence of Q' 7V one finds the additional non-vanishing elements of Af: 



(Af) 17 = -(Af) 47 = -iV, (Af 



g- ' , V ,27 V JO, g, 

The last result can be traced back to the non-vanishing anomalous dimension 7j of the 
weak current at the two-loop level in the HV scheme. Consequently in this scheme the 
operators Q-jy and Q-ja carry anomalous dimensions at the two-loop level. Similarly the 
anomalous dimensions of the operators Q\—Qq receive additional large contributions 
which are related to 7j 7^ in the HV scheme. This feature introduces an unnecessarily 
large scheme dependence at the lower end of the renormalization group evolution. 
In the spirit of []T2] we perform therefore an additional finite renormalization which 
removes these additional contributions to the anomalous dimensions of Qi-Q§ and 
makes the anomalous dimensions of Qjv and Q-ja in the HV scheme to be zero. This 
transformation is equivalent to the following changes in Af of ( |A.2| ): 



.(Af; 



37 



(Af; 



77 



2Ci 



(A.2) 



(Af) 6x6 
(Af) 77 



(Af) 6x6 - 4CVi 6x6 , 
(Af) 77 - 2C F = 0. 



(A.3) 
(A.4) 



The shift in the element (7,7) differs from the shift in (|A.3|) . Therefore, according 
to eq. ( |2.28| ) this finite renormalization modifies also the mixing of the operators Q\- 
Qq with Q' 7V . Using (|2.28|) , the matrix 7W calculated in the NDR scheme and Af 
discussed above we find 7^ in the HV scheme modified by the finite transformations 



( |A.3|) and ( |A.4| ). Here we give only the seventh column of 7#y : 

, ...-m / (In 020 



t8> 



■>£> 



-v (1) 
757 



7g» 



HV 



HV 



HV 



HV 



HV 



HV 



1 

N 




N 



(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 
(A.10) 
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We observe that 7^7 , 7^ and 777'' are the same in both schemes. 
Finally we give the initial conditions in the HV scheme: 

7a s (M w ) , Ajf 7a s (M w ) 



d(M, 



w 



An 



C 2 (M, 



w 1 



1 



6 4vr 



(A.11; 



The other coefficients are given in eqs. (|2.19|) - (|2.24 ) with the constant pieces in ( 2.22 ) 
removed. 



B Evanescent Differences 



As described in section |3.2| , the two-loop anomalous dimension matrix entries 717* - 7^7 
have been extracted from the calculations presented in refs. |S], |TI| [TT 



m eqs. 

In this appendix, we would like to show that all the 7^ - 7g 7 ' ) found here agree with 
the mixings found independently by one of us in ref. p9| , after the differences in the 
conventions used are taken into account. 

The operators responsible for the b — ► se + e~ decay considered in ref. [f29| differ 
from the ones in eq. ( |2.1|) of the present paper only by changes in flavours and by 
trivial normalization factors. The changes in flavours are irrelevant for the values of 
the mixings, as long as the numbers of active up- and down-quarks are kept as free 
variables, which was the case in both papers. In the present paper, we have followed 
the normalization conventions of refs. || |9[ [T0|, 11]. In consequence, the operators 
Qi-Qe in eq. ( |2.1| ) here are four times larger than 0\ - Oq in eq. (3) of ref. pSfl . The 
operator Q' 7V here is two times larger than Og defined in eqs. (3) and (6) of ref. 
Comparing in addition eqs. (|2.10| ) here and (20) in ref. [29| one verifies that 
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7^ 

H7 
H7 



A(B 



1- 
1. 



6. 
.6 



(b.i; 



where the quantities on the l.h.s are the ones from eqs. (|3.1|) -( |3~3"1) and ( p.5|) -( ^9|) of the 
present paper, while B\ z and Bf 3 are taken from eq. (21) of ref. [P| with the following 
substitutions made: 



Q = uQ u + dQ a 




f = u + d 



and £ = 0. (B.2) 



Finally, ABf 3 is the shift required in order to take into account the difference between 
the renormalization schemes used. In the following, we will check that its value is 
exactly what is required for eq. ( |B.1| ) to hold. 

In the present calculation we have used the results of refs. |9|, 0, |TT| obtained within 
the MS scheme and NDR regularization. This was also the scheme used in ref. ||29|1 . 
However, in order to fully specify a scheme, one also has to fix certain convention- 
dependent constants in the definitions of the evanescent operators, whenever such op- 
erators enter the game. Generalizing eq. (15) of ref. |29[, one can write the evanescent 



operators relevant in our calculation in the following form: 



O 



■t+10 



('Y fl ^%rfu'yp])±[i + bi(<i-D)]o i , i 



(B.3) 
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with "— " for i = 1, ...,4 and "+" for i = 5,6. D denotes the dimensionality of 
spacetime. The numbers 6j are arbitrary. Each choice of them corresponds to a different 
renormalization scheme. 
In ref. [E9f all the b^s were set to zero. On the other hand, in the papers [0, 0, IT - 



11 



they differed from zero, and they were fixed by requiring that certain projections of the 
evanescent operators vanish. For i — 1, . . . ,4, the relevant identity was (see eq. (4.3) 
of ref. IDQ 



= E a ^ s (7r(l + 75))^ (7 r (l + 7s)k = 2Tr{| 7 [M7,7p]7r7 [M 7"7 p] 7 T (l + 7b) 
-[1 + 6,(4 - £>)] 7m7t7 /V(1 +75)} = 64(4 - D)(b t - §). 

Eq. (4.8) of ref. |10| gives the corresponding requirement for % = 5,6: 

= ^4(1- 75)^(1 + 7 6 ) fa = 2Tr{| 7[M7l , 7p]7 Vr ] (l+75) + 

+[l + 6 l -(4- J D)] 7M r(l + 75 )} = 32(4 -£>)(&, + §). 

Consequently, the conventions of refs. R |T^, |TTJ correspond to taking 

5 



(B.4) 



(B.5) 
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1 

6' 



6 fi = — 



6' 



(B.6) 



to find the explicit values of 



In the following, we will use the notation of ref. 
AB? 3 for i = 1,...,6. 

According to eq. (16) of ref | 29fl , the contribution of the evanescent operators to 
the two-loop anomalous dimension matrix for the "usual" operators depends on the 
divergences they renormalize (matrix a\l) and on their finite one-loop matrix elements 
(matrix a^). The matrix a]l is insensitive to the 0(4 — D) parts of the evanescent 
operators, i.e. it is insensitive to shifts in the numbers b^. However, the finite one-loop 
matrix elements of the evanescent operators depend on these numbers. If one changes 
bi from zero to any nonzero value, then the corresponding change in the one-loop matrix 
element of O^ 10 is proportional to the product of bi and the divergent part of the one- 
loop matrix element of the "normal" operator that is multiplied by bi in the definition 
of Ofl w . Explicitly: 

(B.7) 



lAa (i+10)j 



<H±k). 



Consequently, according to eqs. (8), (16), (20) and (21) of ref. 

A^ 3 = 16tt 2 £ aH k+10) Bl 3 (±b k ), * = 1 

fe=i 



, one gets 
6 



(B.8) 



with "— " for i — 1, ... ,4 and "+" for i = 5,6. Taking the numbers bi from eq. (B.6) 
above, and the matrices a 11 and B l from eqs. (17) and (21) of ref. respectively, 



one obtains the final result for ABf 3 



2 4 8 4 — 4 20 — 

-Nc F Q u , -c f Q u , -c F Q d , -c F Q + -Nc F Q d , 0, —c F Q 



(B.9) 



We neglect irrelevant terms of order C((4 — D) ) in eqs. (tBJjn 



and 



) here. 
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which is exactly what is needed for eq. ( |B.1| ) to hold. 

We have thus shown that the mixings of the four-quark operators into the quark- 
lepton operator found in ref. are in agreement with the ones found here on the 
basis of the results of refs. || [U], , as long as the differences in the conventions are 



properly taken into account. This suggests that the disagreement between refs. |P9d 
and |56| (see app. B of ref. [^9|) should be resolved in favour of ref. [F29J. However, 



we have to point out that the next-to-leading matching conditions and mixings of the 
four-quark operators from ref. have not been used in a fully consistent manner in 
ref. [p9| |. As long as one wants to define the evanescent operators as in ref. |29| (all 



the numbers 6j vanishing), one has to transform accordingly both the next-to-leading 
terms in the matching conditions for the the four-quark operators Q\ and Q2, and the 
whole 6x6 next-to-leading anomalous dimension matrix of the four-quark operators 
that has been calculated in ref. |J. The author of ref. [29] was not aware of the fact 



that the conventions used in ref. || for the evanescent operators were different from 
his own. In consequence, the numerical results of ref. |29| for the b — > s e + e~ decay 
are somewhat larger than they should be. The change is however phenomenologically 
irrelevant since it never exceeds 3% . 

C Calculation of 7^ and 7^7 in the NDR Scheme 



The procedure described below constitutes a synthesis of subsection 5.1 of refs. [10, 11 



The calculation is more complicated than in the pure QCD case of ref. [JTOJ because here 
we only deal with pure 0(a 2 s ) anomalous dimension matrices after a proper rescaling 
of Q 7V - 

In order to solve the problem of closed fermion loops involving 75 in the NDR scheme, 
we have to consider four different bases of operators. All four bases contain the penguin 
operators Q3 - Q' 7V of eq. ( [2.1| ) but differ in the current-current operators Qi and Q 2 . 
The latter are given as follows, 

Basis A: = {s a u/3) v _ A {u/3d a ) v _ A , Q ( 2 u) = {su) v _ A (ud) v _ A . (C.l) 

Basis B: Q? = (s a d p ) v _ A (d p d a ) y _ A , = (sd) v _ A (dd) y _ A . (C.2) 

Basis C: = (sd) v _ A (uu) v _ A , = {s a dp) v _ A (u p u a ) v _ A . (C.3) 

Basis D: Q { fi = (sd) v _ A (dd) y _ A , Q ( 2 d) = (s a dp) v _ A (df>d a ) v _ A ■ (C.4) 



The basis A is the standard basis of eq. ( [2.1|) and the basis B is an auxiliary basis 
needed for the solution of the problem. The bases C and D are simply Fierz conjugates 
of Qi and Q2 in A and B, respectively. Evidently, 

Q[ d) = Q { 2 ] , Q? = Q? ■ (C5) 

Let us next denote by [Qi]i and [Qj] 2 the result of the type-1 (closed fermion loop) 
and type-2 insertions (no closed fermion loop) of an operator Qi, respectively. Then 
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the results for penguin contributions to the row entries of for Q 3 and Q4 can be 
written as follows 



[Q4 P 



u[Q i 2 ) ]i + d[Q i 2 d) ] 1 + 2[Q 



(d.y 



i + 2[Qi a j2, 



1 J 

(dh 



2 • 



(C.6) 
(C.7) 



For the calculation of the entries in the first six columns it is not necessary to introduce 
the bases B and D because in pure QCD no distinction is made between Q^l and Q^l- 

However, here the seventh column 
ITU the bases B and D have to be 



This was in fact the case considered in ref. [flQf| 
involves couplings to the photon and as in ref. 
introduced. 

There are six independent entries in eq. ( |C.6| ) and eq. (|C.7[ ) which have to be found. 
In the following four steps the quoted results for these entries are given for the seventh 
column only. 
Step 1: 

[Q[% and [Q { 2 % can be calculated without any problems by evaluating the diagrams 
of fig. |j. The contributions to the seventh column are given by 



4QA 



C h 



\\Q\ 



(«)i 



-NC F . 



(C.8) 
(C.9) 



Step 2: 

q! u) and Q { 2 ] receive contributions only from type-1 insertions. This allows to find 
[Q\ ]x and [Q^ i ]i by comparing the two-loop anomalous dimension matrices calculated 
in the bases C and A using the relation 



(7 {1) )f = (7 w )^+[Af,7 w j+2^ Af 



where 



Af 



FJc - l r \A 



(CIO) 
(C.ll) 



with f defined in eq. Q2.2S ). 

Because the insertion in the current-current diagrams are identical for these two 
bases, only penguin diagram contributions to and Af enter this relation. On the 
other hand 7^ is the full one-loop matrix discussed in section 3T . A simple calculation 
of finite terms in one-loop penguin diagrams gives 



Af 



1 

3 



/ \ 
1 


v ° y 



/ N \ 
1 


V I 



p 



(C.12) 



where 



P = (0,0,-l/iV,l,-l/iV,l,0), 
P = (0,0,0,0,0,0,1). 



(C.13) 
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In summary, we find 



(«) 



3* 1 



N 2 



Q 



(«) 



*°<1-JV). 
81 V iV 7 



(C.14) 



Step 3: 

The contribution of [Qi ]i to the seventh column can be related to the corresponding 
contribution of [Qi ]i found in step 2. Inspecting the diagrams of fig. ^ we obtain 



Mi 



[Qi 



(«)i 



4 



(iV 2 



(C.15) 



Step 4: 

The calculation of [Q 2 ]i proceeds as follows. Since Q% receives contributions from 
both type-1 and type-2 insertions we can write 



[Q-2 12 • 

with the last entry calculated in step 1. 

In order to find [Q2 ] p we compare 
calculated in the bases D and B using the relation 



(C.16) 



In order to find [Q2 ]p we compare the two-loop anomalous dimension matrices 



(7 



(7 



+ 



(Af) (d) ,(7 (0) ) (d) l +2/3 Af (d) 



(C.17) 



where the index d indicates that now the auxiliary bases D and B are considered. The 
matrix (^t ))^) differs from (7^) only in the first row and the entry (2,7). The first 
row in (■yWy^ equals the second row and 



( 7 gi)W = ( 7 g))W = 5 (J v + l). 



(CIS 



7W and (Af)( d ) involve again only the penguin contributions. From one- loop penguin 
diagrams we find 



(Af) 



(d) 



1 

3 



/ 1 \ 

-1 


I I 



P-- 9 (N-1) 



( 1 \ 

-1 


V J 



p. 



(C.19) 



In this way the last two terms in eq. (|C.17|) can be evaluated. Since both types of 
insertions of Qi have been calculated in steps 1 and 3, the element 1 can finally 

be extracted from eq. ( |C.17| ) when in addition the relation ( |C.5| ) is used. The result is 



Q 



(d) 



5i ( i-iv). 

81 K N ' 



(C.20) 



Putting all this together and using eqs. (|U.6|) and (|C.7|) we find the results for 737 
and 747* for the NDR scheme quoted in eqs. ( |3.6| ) and ( |3.7|) , respectively. 
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D Collection of Numerical Input Parameters 



Quark Masses 

m t = 120 - 250 GeV m b = 4.8 GeV 

QCD and Electroweak Parameters 



m c = 1.4 GeV 





= 200 MeV : 


a s (M z ) 


= 0.109 




= 300 MeV: 


a s (M z ) 


= 0.116 




= 400 MeV: 


a s {M z ) 


= 0.122 


a 


= 1/128 


G F 


= 1.16639- 10~ 5 GeV 


2 &w 


= 0.23 


M w 


= 80.0 GeV 



-2 



CKM Elements and Bk Ranges 



Range (a) : 
Range (b) : 
Range (c) : 



\V US \ 

\Vcb\ 

\v*\ 
\v ch \ 



■ 0.22 

0.040 ± 0.005 
0.040 ± 0.002 
0.043 ± 0.002 



\v ub /v cb \ 

\Vub/Vcb\ 

\v ub /v cb \ 



K Decays and K° - K° Mixing 



7r°e- 



BR(K + 
t Kl = 5.17- 10- 8 s 
r K+ = 1.237- 10- 8 s 
e = 2.258 • 10~ 3 



= 0.0482 



AM k 



0.08 ±0.02 
0.08 ±0.01 
0.09 ±0.01 



498 MeV 
161 MeV 
3.5 • 10~ 15 GeV 



Bk 
Bk 
B K 



V3 



0.70 ±0.20 
0.75 ±0.05 
0.55 ±0.05 



1.1 

0.57 

0.36 



B° - B° Mixing 



0.72 ±0.08 



r B^ } F B 



200 ± 30 MeV 



i] 2B = 0.55 
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E Tables 



Table 1: Wilson coefficients at [i — 1 GeV for m t = 150 GeV. In our approach 
y 1 = y 2 = holds. 





A m = 0.2 GeV 


A _ = .3GeV 


= 0.4 GeV 


Scheme 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Zl 


-0.583 


-0.394 


-0.474 


-0.710 


-0.483 


-0.602 


-0.848 


-0.584 


-0.769 


Z2 


1.310 


1.196 


1.248 


1.398 


1.254 


1.337 


1.500 


1.325 


1.462 


Z3 


0.004 


0.008 


0.004 


0.005 


0.013 


0.008 


0.007 


0.021 


0.015 


Zi 


-0.010 


-0.023 


-0.011 


-0.013 


-0.034 


-0.018 


-0.018 


-0.054 


-0.029 


Zb 


0.003 


0.006 


0.003 


0.004 


0.007 


0.004 


0.006 


0.009 


0.006 


z 6 


-0.010 


-0.023 


-0.010 


-0.015 


-0.035 


-0.016 


-0.020 


-0.055 


-0.026 


z 7V /a 


-0.015 


-0.008 


0.009 


-0.026 


-0.037 


0.000 


-0.037 


-0.071 


-0.009 



2/3 


0.027 


0.022 


0.024 


0.033 


0.028 


0.032 


0.041 


0.035 


0.042 


2/4 


-0.047 


-0.043 


-0.045 


-0.055 


-0.051 


-0.054 


-0.063 


-0.059 


-0.064 


2/5 


0.011 


0.004 


0.012 


0.012 


0.001 


0.015 


0.013 


-0.008 


0.018 


2/6 


-0.078 


-0.071 


-0.065 


-0.101 


-0.098 


-0.086 


-0.129 


-0.141 


-0.117 




0.554 


0.718 


0.711 


0.524 


0.708 


0.700 


0.498 


0.698 


0.688 


Vta/oi 


-0.579 


-0.579 


-0.579 


-0.579 


-0.579 


-0.579 


-0.579 


-0.579 


-0.579 



Table 2: Wilson coefficient z-jyja for m t = 150 GeV for various values of ji. 





A_ = .2GeV 


A— = 0.3 GeV 


A— = 0.4 GeV 


^[GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


0.8 


-0.030 


-0.020 


0.010 


-0.052 


-0.067 


-0.007 


-0.075 


-0.127 


-0.021 


0.9 


-0.021 


-0.013 


0.010 


-0.037 


-0.049 


-0.002 


-0.053 


-0.094 


-0.015 


1.0 


-0.015 


-0.008 


0.009 


-0.026 


-0.037 


0.000 


-0.037 


-0.071 


-0.009 


1.1 


-0.009 


-0.005 


0.007 


-0.017 


-0.028 


0.002 


-0.025 


-0.055 


-0.005 


1.2 


-0.005 


-0.004 


0.005 


-0.010 


-0.023 


0.002 


-0.015 


-0.044 


-0.002 
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Table 3: Wilson coefficient yrv/oi for m t = 150 GeV for various values of \x. 





A _ = .2GeV 


= 0.3 GeV 


Aj^s = 0.4 GeV 


//[GeV] 


LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


0.8 


0.558 


0.722 


0.715 


0.529 


0.712 


0.704 


0.503 


0.700 


0.691 


0.9 


0.556 


0.720 


0.713 


0.526 


0.710 


0.702 


0.500 


0.700 


0.690 


1.0 


0.554 


0.718 


0.711 


0.524 


0.708 


0.700 


0.498 


0.698 


0.688 


1.1 


0.552 


0.717 


0.709 


0.523 


0.706 


0.698 


0.496 


0.696 


0.686 


1.2 


0.551 


0.715 


0.708 


0.521 


0.705 


0.696 


0.494 


0.694 


0.684 



Table 4: Wilson coefficients y 7 y ja (at /i — 1 GeV) and y-jAjot for various values of m t . 





y 7V /a 


yiAjoi 








0.2 GeV 




0.3 GeV 




0.4 GeV 






GeV] 


No QCD 


LO 


NDR 


LO 


NDR 


LO 


NDR 




130 


0.763 


0.522 


0.686 


0.493 


0.676 


0.466 


0.666 


-0.464 


150 


0.795 


0.554 


0.718 


0.524 


0.708 


0.498 


0.698 


-0.579 


170 


0.823 


0.583 


0.747 


0.553 


0.737 


0.527 


0.727 


-0.703 


190 


0.849 


0.609 


0.774 


0.580 


0.763 


0.553 


0.753 


-0.836 


210 


0.874 


0.634 


0.799 


0.605 


0.788 


0.578 


0.778 


-0.978 


230 


0.898 


0.658 


0.823 


0.629 


0.812 


0.602 


0.802 


-1.129 


250 


0.921 


0.681 


0.846 


0.652 


0.836 


0.626 


0.825 


-1.290 



Table 5: PBE coefficient P of y 7V for various values of Aj^g and /i. In the absence of 
QCD P = 8/9 sm 2 9 w \n(M w /m c ) = 0.827 holds universally. 





Po 


AMsfGeV] 


//[GeV] 


LO 


NDR 


HV 




0.8 


0.487 


0.725 


0.714 


0.2 


1.0 


0.482 


0.719 


0.709 




1.2 


0.477 


0.715 


0.704 




0.8 


0.446 


0.711 


0.699 


0.3 


1.0 


0.440 


0.705 


0.693 




1.2 


0.434 


0.700 


0.687 




0.8 


0.408 


0.694 


0.681 


0.4 


1.0 


0.401 


0.690 


0.676 




1.2 


0.396 


0.684 


0.670 
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Table 6: BR(K L -> 7r°e + e -) dir • 10 12 at n = 1.0 for V cb = 0.041, V ub /V cb = 0.10, 
B K = 0.7 and various values of A^g and m t . The labels (I) and (II) denote the two 
different solutions for the CKM phase 5. 





LO 


NDR 


HV 


Ams [ GeV] 


m t [GeV] 


(I) 


(II) 


(I) 


(II) 


(I) 


(II) 




140 


4.01 


2.95 


5.48 


4.02 


5.40 


3.97 




160 


4.70 


2.44 


6.08 


3.15 


6.01 


3.11 


0.2 


180 


5.18 


2.20 


6.41 


2.72 


6.35 


2.70 




200 


5.61 


2.07 


6.69 


2.46 


6.63 


2.44 




220 


6.02 


1.98 


6.96 


2.29 


6.91 


2.28 




240 


6.41 


1.93 


7.24 


2.18 


7.20 


2.17 




140 


3.79 


2.79 


5.38 


3.95 


5.29 


3.89 




160 


4.49 


2.33 


5.98 


3.10 


5.91 


3.06 


0.3 


180 


5.00 


2.12 


6.32 


2.69 


6.25 


2.66 




200 


5.45 


2.01 


6.61 


2.44 


6.55 


2.41 




220 


5.87 


1.93 


6.90 


2.27 


6.84 


2.26 




240 


6.28 


1.89 


7.18 


2.16 


7.13 


2.15 




140 


3.61 


2.65 


5.27 


3.87 


5.18 


3.81 




160 


4.32 


2.24 


5.89 


3.05 


5.80 


3.01 


0.4 


180 


4.84 


2.06 


6.24 


2.65 


6.16 


2.62 




200 


5.30 


1.95 


6.54 


2.41 


6.47 


2.38 




220 


5.74 


1.89 


6.83 


2.25 


6.77 


2.23 




240 


6.17 


1.86 


7.12 


2.15 


7.07 


2.13 
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F Feynman Diagrams 



Figure 1: One- loop current- current and penguin diagrams in the full theory. 



Figure 2: The three basic ways of inserting a given operator into a four-point function: 
(a) cur rent- current-, (b) type-1 penguin-, (c) type-2 penguin- insertion. The 4- vertices 
"® <g>" denote standard operator insertions. 



Figure 3: One-loop current- current and penguin diagrams contributing to 7*- - 1 in the 
effective theory. The unlabeled wavy lines denote either a gluon or a photon. The 
meaning of vertices is the same as in fig [| Possible left-right or up-down reflected 
diagrams are not shown. 
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(1) 



,1 = 1, 



, 6. Square ver- 
c), respectively. 



Figure 4: Two-loop penguin diagrams contributing to x i7 
tices stand for type-1 and type-2 penguin insertions as of figs. §(b) an d 
Possible left-right reflected diagrams are not shown. The numbering of the diagrams 
corresponds to the notation in ref. JlO 
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G Figures 



Figure 5: Wilson coefficients \y-jyja^ and |?/7A/a| 2 as a function of m t for Aj^t 
0.3 GeV at scale p = 1.0 GeVl. 



Figure 6: Allowed ranges for BR{K L — > 7r°e + e _ ) dir as a function of m t for A^ = 
0.3 GeV at scale \i = 1.0 GeV. The labels (I) and (II) refer to the two possible solutions 
for the CKM phase 5. Parameters \V c b\, \V u b/V c b\ and Bk were varied within present 
experimental limits as given in ( |5.12| ). The light and dark shaded areas correspond 
to the LO and NLO results, respectively. The solid line inside the dark shaded area 
describes the additional restriction from B°-B mixing (see the text). 



Figure 7: The same as in fig. |6|, but for the parameter range given in (|5.13|) . 
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Figure 8: The same as in fig. |6|, but for the parameter range given in ( 5.14 ). 



Figure 9: BR(K L — > 7r°e + e )di r /(ImA t ) 2 as a function of m t for various values of 
Af^g at scale /i = 1.0 GeV. 



32 



References 

[1] B. Winstein and L. Wolfenstein, Rev. Mod. Phys. 65 (1993) 1113. 

[2] J. L. Ritchie and S. G. Wojcicki, Rev. Mod. Phys. 65 (1993) 1149. 

[3] L. Littenberg and G. Valencia, Ann. Rev. Nucl. Part. Phys. 43 (1993) 729. 

[4] M. Kobayshi and K. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 

[5] A. Wagner (NA31), talk presented at the European Physics Symposium on High Energy 
Physics, Marseille, July 1993 . 

[6] L. K. Gibbons et al, Phys. Rev. Lett. 70 (1993) 1203. 

[7] G. Altarelli, G. Curci, G. Martinelli, and S. Petrarca, Nucl. Phys. B187 
(1981) 461. 

[8] A. J. Buras and P. H. Weisz, Nucl. Phys. B333 (1990) 66. 

[9] A. J. Buras, M. Jamin, M. E. Lautenbacher, and P. H. Weisz, Nucl. Phys. B370 

(1992) 69; addendum ibid. Nucl. Phys. B375 (1992) 501. 

[10] A. J. Buras, M. Jamin, M. E. Lautenbacher, and P. H. Weisz, Nucl. Phys. B400 

(1993) 37. 

[11] A. J. Buras, M. Jamin, and M. E. Lautenbacher, Nucl. Phys. B400 (1993) 75. 

[12] A. J. Buras, M. Jamin, and M. E. Lautenbacher, Nucl. Phys. B408 (1993) 209. 

[13] M. Ciuchini, E. Franco, G. Martinelli, and L. Reina, Phys. Lett. 301B (1993) 
263. 

[14] M. Ciuchini, E. Franco, G. Martinelli, and L. Reina, Ecole Normale Superieure 
preprint LPTENS 93/11. 

[15] F. J. Gilman and M. B. Wise, Phys. Rev. D20 (1979) 2392. 

[16] F. J. Gilman and M. B. Wise, Phys. Rev. D21 (1980) 3150. 

[17] C. O. Dib, I. Dunietz, and F. J. Gilman, Phys. Lett. 218B (1989) 487. 

[18] C. O. Dib, I. Dunietz, and F. J. Gilman, Phys. Rev. D39 (1989) 2639. 

[19] J. M. Flynn and L. Randall, Nucl. Phys. B326 (1989) 31. 

[20] A. G. Cohen, G. Ecker, and A. Pich, Phys. Lett. 304B (1993) 347. 



[21] A. Pich, CERN preprint CERN-Th-7114-93; Bulletin Board: |hep-ph 9312297 . 

[22] M. Misiak and Z. Pluciennik, Phys. Lett. 268B (1913) 429. 

[23] G. Buchalla, A. J. Buras, and M. K. Harlander, Nucl. Phys. B349 (1991) 1. 

[24] M. K. Gaillard and B. W. Lee, Phys. Rev. Lett. 33B (1974) 108. 

[25] G. Altarelli and L. Maiani, Phys. Lett. 52B (1974) 351. 



33 



[26] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, JEPT 45 (1977) 670. 

[27] B. Guberina and R. D. Pecchei, Nucl. Phys. B163 (1980) 289. 

[28] J. O. Eeg and I. Picek, Phys. Lett. 214B (1988) 651. 

[29] M. Misiak, Nucl. Phys. B393 (1993) 23. 

[30] G. Buchalla and A. J. Buras, Nucl. Phys. B398 (1993) 285. 

[31] G. Buchalla and A. J. Buras, Nucl. Phys. B400 (1993) 225. 

[32] G. Buchalla and A. J. Buras, Nucl. Phys. B412 (1994) 106. 

[33] S. Herrlich and U. Nierste, Technical University Munich preprint TUM-T31- 



49/93; Bulletin Board: hep-ph 9310311 



[34] A. J. Buras, M. Jamin, and P. H. Weisz, Nucl. Phys. B347 (1990) 491. 

[35] W. A. Kaufman, H. Steger, and Y. P. Yao, Mod. Phys. Lett. A3 (1988) 1479. 

[36] G. Buchalla, A. J. Buras, and M. K. Harlander, Nucl. Phys. B337 (1990) 313. 

[37] A. Datta, J. Frohlich, and E. A. Paschos, Zeitschr. f. Physik C46 (1990) 63. 

[38] J. M. Flynn, Mod. Phys. Lett. A5 (1990) 877. 

[39] G. W. Kilcup, S. R. Sharpe, and R. Gupta, University of Washington preprint 
UW/PT-93-25. 

[40] W. Bardeen, A. J. Buras, and J.-M. Gerard, Phys. Lett. 211B (1988) 343. 

[41] J. Prades et al, Zeitschr. f. Physik C51 (1991) 287. 

[42] A. S. Kronfeld and P. B. Mackkenzie, Ann. Rev. Nucl. Part. Sci. 43 (1993) 30. 

[43] D. G. Cassel, talk given at Max-Planck Institute, Munich, Dec. 15, 1993 . 

[44] D. A. Harris et al, Phys. Rev. Lett. 71 (1993) 3918. 

[45] K. E. Ohl et al, Phys. Rev. Lett. 64 (1990) 2755. 

[46] G. Cella, G. Ricciardi, and A. Vicere, Phys. Lett. 258B (1991) 212. 



34 



12.0 



„ 10.0 

CM 
i 

o 

— 8.0 



o 

T 4.0 




0C 2.0 
CD 



0.0 



120 140 160 180 200 220 240 

m t [GeV] 




